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development  of  an  operational  warm  fog  dispersal  system  using  momentum- 
driven  ground-based  heat  sources.  In  previous  investigations,  the  wind  was 
in  the  same  direction  as  the  jet  (coflowing).  The  present  investigation  is 
concerned  with  a jet  opposite  in  direction  to  the  ambient  wind  (counterflowing 
or  opposing  jet). 

v.A  model  for  calculating  the  dynamic  characteristics  of  a cold  counter- 
flowing  jet  has  been  extended  to  take  account  of  the  effect  of  heating  upon  the 
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tend  to  be  greater  than  the  experimental  results  for  the  heated  jet  at  large  jet 
speeds  relative  to  ambient  velocity. 

The  dependence  of  jet  length,  axial  velocity,  and  axial  temperature  upon 
initial  jet  temperature  is  obtained  for  both  planar  and  round  jets,  and  the 
results  are  consistent  with  those  obtained  for  a submerged  nonisothermal  jet. 
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Interaction  of  a Turbulent 
Planar  Heated  Jet  with  a 
Counterflowing  Wind 


1.  INTROIM  ( TION 


To  aid  in  the  development  of  an  operational  warm  fog  dispersal  system  (VVFDS), 
experimental  and  theoretical  studies  have  been  made  of  the  characteristics  of 
ground-based  heated  jets  for  various  combinations  of  heat  and  thrust  under  differ- 
ent wind  conditions.  The  dynamic  characteristics  and  trajectories  of  buoyant 

12  3 

heated  jets  in  a coflowing  wind  have  been  presented  in  previous  reports.  ’ ’ In 
references  1 and  2,  the  characteristics  of  planar  and  round  jets  were  obtained;  in 
reference  3,  a method  was  developed  for  taking  into  account  the  attachment  of  the 
jet  to  the  ground  for  a considerable  portion  of  the  jet  trajectory  (ground  effect). 

The  dynamic  characteristics  of  non-buoyant  coflow'ing  jets  are  well  known  and 

4 

were  presented  in  detail  by  Abramovich.  The  principle  problem,  therefore,  in 
determining  the  trajectory  of  a heated  coflowung  jet  is  to  properly  evaluate  the 

(Received  for  publication  2fi  September  1977) 

1.  Klein,  M.M.,  and  Kunkel,  B.A.  (1975)  Interaction  of  a Buoyant  Turbulent 

Planar  Jet  With  a Co-Flowing  Wind,  A FC RL-TR-75-()368. 

2.  Klein.  M.M.,  and  Kunkel,  B.A.  (1975)  Interaction  of  a Buoyant  Turbulent 

Round  Jet  With  a Co-Flowing  Wind.  AFCRI.-TR-75-0581. 

3.  Klein,  M.M.  (1977)  A Method  for  Determining  the  Point  of  Lift-Off  and 

Modified  Trajectory  of  a Ground-Based  Heated  Turbulent  Planar  Jet  in  a 
Co-Flowing  Wind,  A FC.I.-TR -77-0033. 

4.  Abramovich,  G.  N.  (1963)  The  Theory  of  Turbulent  Jets,  The  MIT  Press, 

Cambridge,  Mass. , Chaps.  4 and  9. 


7 


buoyant  force.  The  situation  for  the  counterflowing  jet,  as  presented  by 
4 

Abramovich,  is  considerably  less  satisfactory  for  both  theory  and  experiment 
because  the  calculations  that  must  account  for  the  regions  of  counterflow  are 
quite  specialized  and  not  easily  adaptable  to  the  buoyant  jet.  Recently,  however, 
a simplified  model  was  presented  by  Sekundov0  for  a counterflowing  round  jet,  and 
it  can  be  extended  in  a straightforward  w»ay  to  take  account  of  buoyancy  and  obtain 
the  jet  trajectories.  An  important  feature  of  the  Sekundov  model  is  the  use  of  a 
finite  wall  to  help  simplify  the  equations  of  motion.  The  results  for  an  open  jet 
are  then  obtained  by  making  the  wall  diameter  arbitrarily  large. 

Because  Sekundov  limits  himself  to  the  unheated  incompressible  jet,  the 
model  must  first  be  extended  to  take  account  of  density  variation.  In  addition,  an 
energy  equation  must  be  developed  to  help  determine  the  temperature  distribution 
along  the  jet  axis.  The  effect  of  buoyancy  upon  the  jet  characteristics  will  be 
presented  in  a subsequent  report.  Since  the  jets  in  a warm  fog  dispersal  system 
merge  a short  distance  downstream  of  the  jet  nozzles,  the  analysis  has  been 
confined  to  the  planar  case. 


2.  JIT  t.KOMKTR) 

A schematic  sketch  of  the  flow  pattern  for  a planar  jet  in  a counterflowing 
wind  is  shown  in  Figure  1.  The  jet  issues  with  a velocity  u and  starts  mixing 
with  the  counterflowing  stream  of  velocity  ua,  resulting  in  a boundary  layer  that 
grows  inward  toward  the  center  of  the  jet.  This  region,  called  the  initial  section, 
terminates  when  the  boundary  layer  reaches  the  center  of  the  jet.  The  surface 
yjix),  where  x is  horizontal  position  and  y the  vertical  coordinate,  separates  the 
regions  of  unperturbed  and  perturbed  flow  in  this  initial  section. 

The  velocity  profile  remains  constant  in  the  next  interval,  known  as  the  main 
section.  The  surface  y.,(x)  separates  the  region  in  which  the  flow  has  retained  its 
original  direction  from  that  in  which  the  flow  has  been  reversed.  Thus,  at  y the 
velocity  is  zero.  Analogous  to  the  surface  y ^ in  the  initial  section,  the  surface 
y-5(x)  separates  the  regions  of  perturbed  and  unperturbed  flow'  for  the  opposing 
stream.  The  line  y3  is,  therefore,  a streamline  of  the  flow. 

Because  of  the  constricting  effect  of  the  walls,  the  unperturbed  velocity  uc  in 
the  main  region  is  slightly  greater  than  u^.  Experiment  indicates  that  the  flow 


5.  Sekundov,  A.N.  (1969)  The  propagation  of  a turbulent  jet  in  an  opposing 

stream,  in  Turbulent  Jets  of  Air,  Plasma,  and  Real  Gas,  Consultants 
Bureau,  New  York,  A Division  of  Plenum  Publishing  Corp.  , N.  Y.  10011. 

6.  Kunkel,  B.A.  (1975)  Heat  and  Thrust  Requirements  of  a Thermal  Fog 

Dispersal  System.  AFCRI.-TR-75-0472. 
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Figure  1.  Schematic  Representation  of  Geometry  of  Counterflowing  Jet  and 
Velocity  and  Temperature  Profiles 

regime  in  a counterflowing  jet  can  be  separated  into  two  zones.  In  zone  I,  consist- 
ing of  the  initial  and  main  sections,  the  pressure  is  fairly  constant  and  the  surface 
y.,  grow's  at  a constant  rate.  In  zone  II,  the  pressure  starts  to  rise  rapidly,  the 
surface  y.,  starts  to  decrease,  and  the  final  retardation  of  the  jet  flow  is  attained. 
.Since  at  this  point  the  axial  velocity  u is  comparable  in  magnitude  to  the  stream 

velocitv,  we  mav,  for  convenience,  here  take  u as  numerically  equal  to  u . 

m J M a 

Because  of  the  large  pressure  rise  in  zone  II,  the  velocity  and  .temperature 
distributions  in  this  regime  cannot  be  accurately  calculated.  How'ever,  since  these 
quantities  do  not  differ  greatly  from  their  free-stream  values,  simple  interpolation 
may  be  used  where  required. 

:i.  HI. mill  HI)  TKMI’Kll  \ II  HI  I’KOKII.KS 

Following  Sekundov  we  shall  utilize  a linear  profile  for  the  transverse  velocity 
distribution  and  will  assume,  in  addition,  that  the  temperature  profile  is  a linear 
function  of  the  velocity  profile.  The  use  of  linear  profiles  considerably  simplifies 
the  analysis,  and  should,  as  is  usually  the  case  with  integral  methods,  have  only 
a weak  effect  upon  the  final  results.  For  the  region  between  the  axis  and  y„,  we 
write  the  velocity  and  temperature  distributions  in  the  form 


u . 

1 - V;  Vr 

u *2 

m 


(1) 


u *•  u 


u 


m 


c_ 

n 

c 


(2) 


where  subscript  m refers  to  values  along  the  axis  and  subscript  c to  the  unper- 
turbed flow  above  v.^ . 

At  v y9,  u 0 so  that  the  temperature  T.,  at  y,;  is  given  by 


(2) 


For  the  region  between  y9  and  y., , the  profiles  have  the  form 


(4) 


(5) 


i.  \\i\i  nisimm  non  m ui.oun  win i:\iim  h \ ri  m 


hi  Inilial  Section  \n<il  \ -»i  for  Mass  ,tml  Momentum 

Conservation  of  mass  and  momentum  between  the  beginning  of  the  initial 
section,  O,  and  section  V to  the  right  of  zone  II  yields 
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wh  re  is  the  tensity,  p the  pressure,  and  the  subscript  o refers  to  initial  jet 

values,  since  y , is  a streamline  of  the  flow,  the  pressure,  density,  and  velocity- 
in  the  region  above  y^  are  related  by  the  Bernoulli  equation  for  channel  flow, 

f ^ + 4-  const  . (8) 

In  the  development  of  the  conservation  equations,  we  will  be  interested  in  the 
limiting  form  of  these  equations  as  the  wall  width  becomes  large.  In  the  region 
above  v.,,  a typical  quantity,  that  is,  density-  pc,  may  be  replaced  by  pa,  except 
where  the  limiting  form  (pa  - pc)b  occurs.  We  now  show  that  the  only  limiting 
form  of  consequence  is  that  containing  the  velocity,  the  others  being  negligible. 
To  evaluate  the  integral  in  Kq.  (8),  we  utilize  the  adiabatic  law 

-5-  - const  (9) 

. “V 


where  y is  the  ratio  of  specific  heats,  yielding 


o •> 


(10) 


Since  ,,  is  close  to  p , we  mr.v  express  in  the  form 
r c a ■ *c 

Pc  Pa  - (pa  " pc' 

and  expand  p to  fir^t  lrder  in  the  small  difference  p - p(  to  yield 


a , 2 2.  1 a 

— (u  - u 1 ' — 
pens  p 
Ka  a 


( u 


u ) 
a 


From  the  adiabatic  law,  Fq.  (P),  the  pressure  difference  is  given  by 


(11) 


(12) 


while  from  the  perfect  gas  lav. 


the  temperature  difference  is 


(7 


(14) 


For  a typical  case,  u 20  m/sec  (moderately  high  speed  jet),  p 1 atm, 

SL  o 3 

T&  287°K,  0.00123  gm/cm'  , and  the  coefficient  on  the  right  side  of  Eq.  (11) 

has  the  value 

1 P o 9 

± — u = 0. 00348  . 

y p a 

1 a 

Thus,  the  density  difference  and,  from  Eqs.  (12)  and  (14),  the  pressure  and 
temperature  difference  are  all  small  compared  to  the  velocity  difference.  We 
may  therefore  replace  p and  T by  p and  T and  write  the  pressure  difference 
in  the  form 


P - P 
ra  c 


(u“ 


u2) 

a 


(15) 


Because  of  the  virtual  constancy  of  the  pressure  in  zone  I,  the  integrals  in 
Eqs.  (fi)  and  (7)  can,  by  use  of  the  perfect  gas  law,  be  easily  evaluated.  However, 
the  logarithmic  forms  introduced  by  the  integrals  makes  the  equations  awkward  to 
handle  in  subsequent  analysis.  Since  we  are  interested  in  small  to  moderate  tem- 
perature and  density  differences,  we  will,  for  convenience  and  simplification, 
utilize  an  arithmetic  mean  for  the  density  difference  and  write 


average 


p 1 (p2  + Pm’ 

tZ  al  2 — o 


0 < y < vr 


average  - — a 0 

p o 2 


1 (p2  4 P a) 


y0  £ y s y„ 


The  mass  and  momentum,  Eqs.  (6)  and  (7),  now  take  the  form 


(16) 

(17) 
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At  this  point  we  take  all  quantities  as  dimensionless,  referring  to  ambient 

values  for  density,  velocity,  and  temperature  and  to  the  initial  iet  width,  v for 

' o' 

length.  The  mass  and  momentum  equations  may  now  be  written  as 


(y,n  “ y9n> 

(uc  - l)b  p0m  - — ■ , - ao  + y30 


pom  * 3 (y30  * •V20)  + y30  ' 2 <Uc  ' 1>b 


(20) 


(21) 


where  m is  used  for  the  dimensionless  initial  jet  velocity,  u , and  y30  and  y,)f) 
indicate  initial  values  of  y3  and  y9.  Since  y00  1,  Eqs.  (20)  and  (21)  may  be 


20 


solved  for  the  unknowns  v.,,.  and  (u  - l)b  to  vield 

3 0 c 


(i  m p _ 


(m  + 1)  + 1 


(22) 


<Uc  ' 1)b  pom  ' >30  *~T 
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ao  ° ° 
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(23) 


1.2  Main  Section  \nnhsi-.  for  Mass  uml  Momcnluni 

The  conservation  equations  for  the  region  between  a typical  section  I in  the 
main  region  and  section  V are 


y2 


f p u dy  + b - J p u dy  - u^p^tb  - v.,)  0 


y3 


f p u“  dy  - b - J p u2  dy  - pcu“’(b  - yg)  + (p;>  - p^)b 


(24) 


(25) 
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which,  after  use  of  Eq.  (15)  and  integration  of  the  profile  distributions,  may  be 
written  as 


Pmumal  T ' y2(N  - n ~T  ' (uc  - 1)b  + Ny2  ° 


PmUm  X y2  * (N  - l)  X y2  ‘ (uc  ' l)  | + Ny2 


where  N y.^/y,,.  We  may  unite  Eq.  ( 2 H)  in  the  more  compact  form 


[f  * » (‘-if 


8 p u a,  + a 
m m 1 o 


a (u  - l)b 

c 


I'se  of  Eq.  (28)  for  o in  Eq.  (27)  yields 


/ 3 \ 3 i\  H . 

4j-  + N ( 1 - ) N - (N  - 1)  - p u“  -~ 

2 \ 2 / ! ' m m 3 


which  may  he  solved  for  X in  the  form 


a 1 2 a i 

N 2 — — p u f 3—  p u H . 

a mm  a 'mm 


Since  a has  been  obtained  from  its  value  at  the  initial  section,  Eq.  (28),  and 

is  assumed  constant  throughout  the  main  section,  Eqs.  (28)  and  (.82)  provide  the 

distribution  of  axial  velocity  and  the  coordinate  v , if  the  axial  temperature  T is 

2 1 m 

known.  From  Eqs.  (22)  and  (28)  we  may  write  the  velocity  distribution  in  the  more 
explicit  form 


[3(2  - a )m“  + 2 (3  - a )ml  • 1 
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I.  i r<-iii|iiTaliiri'  Viiulvni*  for  Main  Section 


To  obtain  the  temperature  distribution  in  the  main  section,  we  calculate  the 
transport  of  energy  between  sections  1.  and  V and  obtain,  analogous  to  the  trans- 
port of  mass  in  Eq.  (24), 


/ puTdy  -J  puT  dy  - ucpcTc(b  - y3)  -b  . (34) 

° y2 

Expressing  T in  terms  of  temperature  excess  over  ambient,  AT  = T - 1,  and 
making  use  of  the  mass  equation,  we  can  write  Eq.  (34)  in  the  equivalent  form 


/ 


p u AT  dy 


y2 

/ p u AT  dy  = pc  uc(l  - Tc)(b  - y3) 


(3  5) 


The  term  1 - T£  can,  with  the  aid  of  Eqs.  (11)  and  (14),  be  represented  in  terms 

of  the  velocity  difference  Au  = u - 1 bv 

c c J 


CP  Ta  ^ 


(3  (!) 


where  C ^ is  the  specific  heat  of  the  gas,  taken  as  constant.  Utilizing  Eq.  (28), 
expressing  p in  terms  of  T by  the  gas  law,  and  noting  that  the  term  y3  is  negligible 
compared  to  b for  a wide  jet,  the  energy  equation  can  be  written  in  the  form 
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Parallel  to  the  treatment  of  the  density  term  in  Eqs.  (6)  and  (7),  we  replace  1 /T 
by  its  corresponding  arithmetic  mean  p to  give 
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where 
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Performing  the  integration  we  obtain 
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(39) 


(40) 


which,  together  with  Eqs.  (28)  and  (32),  allows  us  to  calculate  temperature, 
velocity,  and  the  coordinate  y3  along  the  main  section.  The  coordinate  y9  occurring 
in  these  equations  is  a simple  linear  function  of  x given  by 


y2  = cx  (41) 

where  c is  a constant  having  the  value  0.  27  for  the  initial  section  and  0.  22  for  the 
main  region. 


Ml. TIIOI)  01  SOM  NON 


Although  we  have  simplified  the  analysis  by  the  use  of  average  values  for 
temperature  and  density,  the  solution  for  the  foregoing  equations  is  quite  complex 
for  the  general  case.  However,  for  the  case  of  small  to  moderate  heating  the 
equations  can  be  further  simplified  and  explicit  solutions  obtained.  The  quantity 
ATm  is  initially  comparable  to  unity  but  decreases  very  rapidly  and  becomes  small 
in  the  earlier  portion  of  the  main  region.  The  term  AT.,  is  always  small  compared 
to  unity.  Thus,  the  coefficient  aQ,  Eq.  (17), 
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(42) 


AT, 


is  close  to  unity.  If  we  write  Eq.  (33)  in  the  form 

P 


f ♦ n (i 


y„  — — [2(1  - a ) f 4m  + 3(1  - a )m  + 3m  1 + 1 (43) 

^2  a o o 


and  use  the  initial  value  of  from  Eq.  (3), 


AT, 
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^To 
1 + m 


(44) 


and  note  that,  for  the  test  data  analyzed,  m is  large  compared  to  unity,  we  can 
write  Eq.  (43)  in  the  approximate  form 

<t>  (m) 

v,  = — ^ (45) 

o 

where 

$o(m)  3m"  + 4m  + 1 + 3(1  + m)  AT^  . (46) 


Utilizing  Eq.  (32)  for  the  value  of  N,  Eq.  (31)  may  be  written  as 
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IT 


+ N 


+ 1 


in  which,  from  Eqs.  (16)  and  (17), 


a,  a is  given  by 
1 o 6 


AT  . 


2 + u 


AT., 


(47) 


(48) 


Since  AT,  AT  /(I  + u ) is  a small  term,  we  can  approximate  Eq.  (48)  by 
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and  reduce  Eqs.  (47)  and  (32)  to  the  form 
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The  velocity  equation  (45)  and  the  energy  equation  (40)  may  now  be  written  as 
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(53) 


The  velocity  equation  (52)  and  the  energy  equation  (53)  may  be  solved  simul- 
taneously in  the  form  or  cubic  equations  for  the  velocity  and  temperature  as  func- 
tions of  jet  position.  A much  simpler  solution  may  be  obtained,  however,  by 
noting  that  AT  drops  off  rapidly  in  the  initial  region  of  the  jet  and,  therefore, 
the  temperature  terms  on  the  right-hand  side  of  Eqs.  (52)  and  (53)  have  only  a 
small  effect  upon  the  velocity  and  temperature  distributions.  We  therefore  neglect 
these  terms  in  obtaining  initial  solutions  and  write  Eqs.  (52)  and  (53)  in  the  form 
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(55) 


whic  h give  .simple  solutions  for  u and  AT  as  functions  of  x.  It  is  convenient  to 
1 m m 

refer  Eqs.  (54)  and  (55)  to  the  end  of  the  initial  section  x„  at  which  u^  m, 

AT  AT  , and  write 

m o 


IB 


(u  + 1)“  (m  + 1) 


2 XH 


(50) 
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(57) 


where  Xjj  is,  by  Eq.  (54),  obtained  from 


<i>0(m) 


°-  27  (m  + l)2  T 


(58) 


Similarly,  the  end  of  the  main  section  x,,  at  which  u 
■"  1 m 


1,  is  given  by 


1 *o(m) 


1 0.22  4T 


(59) 


The  initial  solution  for  the  temperature  obtained  from  Eq.  (57)  may  now  be 
utilized  in  Eqs.  (52)  and  (53)  to  obtain  more  accurate  solutions  for  the  velocity  and 
temperature.  Calculations  indicated  that  the  temperature  correction  has  a small 
effect  upon  the  temperature  distribution  and  a negligible  effect  upon  the  velocity. 

It  may  be  noted  i.hat  Eq.  (57)  for  the  temperature  distribution  is  close  to  tha.  for 
the  coflowing  jet, 


AT  ' u 
m m 


1 


at  small  to  moderate  jet  distances,  but  yields  somewhat  larger  values  in  the  outer 
region  of  the  jet  where  u^  becomes  comparable  to  unity. 
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Utilizing  the  analysis  of  Abramovich  (Ref.  4,  Chap.  9,  pp  391-416),  the  length 
of  zone  II  is  given  in  terms  of  the  length  Xj  of  zone  I by,  approximately, 

X£  - x^  = 0.  176  Xj  . (60) 
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An  analysis  of  the  test  data  obtained  at  Irvine,  California,  during  1974  (Kunkel1’) 
shows  that  the  calculated  jet  length,  x^,  obtained  from  Eq.  (f>0)  is  too  large,  by 
about  a factor  of  2,  compared  to  the  corresponding  experimental  values,  particu- 
larlv  at  the  large  values  of  m.  This  is  in  contradistinction  to  the  results  for  the 
opposing  round  jet  in  which  Sekundov ’ obtains  jet  lengths  that  are  in  good  agree- 
ment with  the  experimental  values.  Part  of  the  discrepancy  may  stem  from  the 
fact  that  we  have  a two-dimensional  flow  that  is  more  sensitive  to  small  changes 
in  velocity  than  an  axisvmmetric  flow.  We  therefore  modify  the  analysis  that  was 
used  previously  to  determine  the  value  of  u(>  from  the  boundary  conditions  at  the 
beginning  of  the  initial  section  and  account  for  the  situation  at  the  end  of  the  main 
section. 

The  continuity  equation  for  the  region  above  y^,  between  the  initial  section  O 
and  an  arbitrary  section  in  the  main  region. 


p u = p U — r- 

Fc  c Fco  co  b 


(b  - y30) 


(fill 


yields  for  large  b 


(u  - l)b  (u  - l)b  + v,  - v,n  (62) 

c co  • 3 ' 30 

in  which  the  density  variation  is  dropped  since  it  is  small  compared  to  the  velocity 
change,  t'se  of  Eqs.  (23)  and  (28)  yields 


The  coordinate  y~  may  be  calculated  from  Eq.  (51)  for  N and  Eq.  (52)  for  y.,  to 
give 
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in  which  we  note  the  term  5AT  (3m  4 1)  is  small  compared  to  unity.  Evaluating 

Eq.  (64)  at  u 1 and  noting  that  here  AT  is  small,  we  get 
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Utilizing  Eq.  (65)  for  y3  and  Eq.  (22)  for  y Eq.  (63)  takes  the  form 


(um  + !)  y* 


1 + AT 


0 j(m) 


AT  , , v T 

'•T1  -nV  " 

m ' 


(60) 


where 
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(67) 


Since  the  right-hand  side  of  Eq.  (66)  is  approximately  1/2  that  of  Eq.  (52), 
the  corrected  value  of  Xj  is  about  1 2 the  value.  This  correction  brings  the  cal- 
culated jet  length  into  closer  agreement  with  the  corresponding  experimental 
value.  However,  the  velocity  u^  is  now  too  small,  especially  in  the  earlier  region 
of  the  jet.  This  situation  is  not  surprising,  since  we  now  have  a solution  that  is 
accurate  at  Xj  but  does  not  satisfy  the  condition  um  m at  x^.  We  therefore 
utilize  a linear  combination  of  the  unmodified  and  the  modified  solutions  that  re- 
duces to  Eq.  (52)  at  xR  and  to  Eq.  (66)  at  x^  that  is. 
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in  which  we  neglect  the  small  term  1/(1  + u^)  in  the  temperature  correction.  It  is 
convenient  to  express  Eq.  (68)  in  terms  of  the  parameters  x^  and  x^  and  write 
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where  xf)  and  Xj  are  related  to  and  <j)  (m)  by 
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(71) 


d^m)  1 + AT  j 


1 4 c T 1 + AT,  /2 

o 1 


and  AT,  is  the  value  of  AT  at  u 1 . 

1 mm 

As  an  alternative  to  the  use  of  Eqs.  (70)  and  (71)  for  x^  and  Xj  in  Eq.  (08), 
we  may  regard  Xjj  and  x^  as  experimentally  determined  parameters  for  calculating 
the  velocity.  Since  the  temperature  distribution  is  given  directly  in  terms  of  veloc- 
ity, Eq.  (53)  may  be  used  unaltered  in  calculating  the  modified  temperature 
distribut  ion. 


7.  KKh'KCT  OK  JKT  UK  (TIM.  I PON  I))  \ Wilt  ( II  \H  \t  11  111'  III  - 01  III 

The  foregoing  analysis  shows  clearly  the  effect  of  heating  upon  the  iet  charac- 
teristics. From  Eq.  (68)  we  see  that  the  velocity  decreases  as  1 1 ~ while 

- 1 2 ° 

Eq.  (55)  indicates  that  AT  decreases  roughly  as  T “.  The  length  of  zone  I,  x., 

® _ i 1 

and,  hence,  the  jet  length  x,,  decrease,  according  to  Eq.  (70),  as  T(  . The  effect 
of  heating  upon  velocity  and  jet  length  for  a round  jet  can  easily  be  deduced  from 
Sekundov's  work  by  noting  that  the  initial  heating  term  introduces  the  density  term 
Po  1/Tq  into  the  term  0((m).  Sekund  ,v's  \ locitv  equation  (3-29)  may  therefore 
be  written  as 
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x“i  (u 


m 


) 
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where  i (u  ) is  a quadratic  in  u 
m _ , , „ ^ m 

roughly  as  T~ 


Thus,  the  velocity  and  iet  length  both  decrease 
Similarly,  it  is  easy  to  show  from  conservation  of  energy  for  a 

round  iet.  analogous  to  Eqs.  (33)  and  (37),  that  the  temperature  drops  off  as 
t-1/2 


As  indicated  by  Sekundov,  1 the  dependence  of  jet  length  upon  temperature  for 
a round  opposing  jet  may  easily  be  deduced  by  noting  that  the  propagation  charac- 
teristics of  zone  I are  similar  to  those  of  a submerged  jet,  and  zone  I is  the  major 
region  of  an  opposing  jet.  In  this  way,  using  conservation  of  momentum  and  energy 
for  a submerged  nonisothermal  jet,  we  obtain. 


for  a round  iet: 


2 2 

xu  ~ p 
m Ko 


AT  u x a 
mm  Ho 


22 


for  a plane  jet; 


OiT  u x - p 
m m o 


from  which  the  previously  obtained  dependence  of  jet  length,  velocity,  and  tem- 
perature upon  initial  jet  temperature  for  both  plane  and  round  jets  is  immediately 
deduced. 
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The  calculated  velocity  and  temperature  distributions  for  several  Irvine  tests 
are  shown  in  Figures  2 and  2 along  with  the  corresponding  experimental  curves. 
These  results  cover  most  of  the  range  of  initial  jet  speed,  initial  jet  temperature, 
and  ambient  wind  encountered  in  the  tests.  In  general,  the  calculated  results  are 
in  fair  to  good  agreement  with  the  experimental  results.  The  calculated  velocity 
drops  off  too  rapidly  in  the  initial  region  of  the  jet,  but  tends  to  approach  the 
experimental  value  at  larger  let  distances.  The  calculated  temperatures  also  drop 
off  somewhat  too  rapidly  in  the  initial  region  of  the  jet,  but  then  they  decrease 
slowly  enough  to  give  temperatures  higher  than  the  experimental  values.  These 
results  are  consistent  with  those  obtained  for  the  coflowing  let  (Klein'). 

The  slower  decrease  of  the  experimental  values  of  temperature  and  velocity 
in  the  early  portion  of  the  jet  are  probably  due  to  incomplete  mixing,  resulting  in 
higher  maximum  values  of  velocity  and  temperature  at  a given  location.  A possible 
factor  contributing  to  the  more  rapid  outer-region  decrease  in  the  experimental 
curves  of  velocity  and  temperature  is  the  buoyant  motion  of  the  jet  that  results  in 
more  rapid  mixing  than  occurs  for  a horizontal  nonbuoyant  ,tet. 

Several  values  of  the  jet  length  x„  obtained  by  the  present  method  are  shown 
in  Figure  4 as  a function  of  m,  along  with  the  corresponding  experimental  results. 

In  addition,  the  experimental  values  obtained  by  Timma'  for  an  incompressible 
unheated  jet  are  presented  to  indicate  the  trend  at  low  values  of  m.  In  order  not  to 
complicate  the  results  due  to  heating,  the  experimental  and  calculated  results  have 
been  restricted  to  initial  excess  temperatures  less  than  30  percent  of  ambient  tem- 
peratures. The  calculated  values  appear  to  follow  the  trend  indicated  bv  the 
results  of  Timma,  while  the  experimental  results  are  somewhat  lower.  The 


7.  Timma,  K.  (1002)  Turbulent  circular  and  planar  lets  in  an  opposing  stream, 
Izv,  Est.  SSK,  Ser . Techn.  i.  I Iz.  -Mat.  Nauk,  11  (No.  4). 
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figure  Comparison  of  Calculated  Temperature  Distributions  Along  Jet 
Axi^  With  Corresponding  Experimental  Curves 


lower  experimental  values  are  probably  due,  in  part,  to  the  buoyancy,  which  tends 
to  decrease  the  jet  length. 
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Figure  4.  Comparison  of  Calculated 
Jet  Length  X2  With  Experimental 
Result  at  Several  Vaiues  of  m 
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The  method  developed  by  Sekundov'  to  calculate  the  dynamic  characteristic  of 
an  incompressible  round  jet  in  an  opposing  stream  has  been  extended  to  the  case 
of  a heated  planar  jet.  The  calculated  velocity  and  temperature  distributions  along 
the  jet  axis  are  in  fair  to  good  agreement  with  the  experimental  results.  In  gen- 
eral, the  calculated  velocity  and  temperature  values  drop  off  too  rapidly  in  the 
early  portion  of  the  jet,  while  the  temperature  decreases  too  slowly  in  the  outer 
region.  The  calculated  jet  lengths  are  consistent  with  experimental  values  for 
unheated  low  velocity  jets,  but  tend  to  be  higher  than  the  experimental  results  for 
heated  high  velocity  jets. 

The  analysis  shows  that  for  a planar  jet  the  velocity  and  temperature  decrease 

_l/2  -1 

as  T ' , while  the  jet  length  drops  off  as  T . For  a round  heated  jet  the  velocity, 

1 /2 

jet  length,  and  temperature  all  decrease  as  T^  . These  results  are  consistent 
with  those  obtained  by  simple  dimensional  analysis  of  the  conservation  laws  for 
momentum  and  energy  for  a submerged  nonisothermal  jet. 
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List  of  Sy  m bols 

b distance  from  outside  wall  to  jet  axis 

c jet  thickness  coefl'icient 

H end  of  initial  section 

m wind  speed  parameter  (u  u ) 

o a 

L,  arbitrary  position  in  main  section 

O initial  position  in  main  section 

p pressure 

T jet  temperature 

T ambient  temperature 

T temperature  on  axis 

AT  temperature  excess  over  ambient  (T  - T ) 
u jet  velocity 

V arbitrary  section  in  main  stream 

x horizontal  position  along  jet 

y vertical  position 

v jet  half-width 

yj  zero  velocity  surface  in  initial  section 
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zero  velocity  surface  in  main  section 

surface  separating  perturbed  and  unperturbed  flows 

gas  density 


y2 

>3 

P 

Subscripts 

a ambient 

c region  above  y^  surface 

m on  jet  axis 

o initial  jet  position 


